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ABSTRACT 

The spatial distribution of persistent spins at zero-temperature in the pure two-dimensional 
Ising model is investigated numerically. A persistence correlation length, £(£) ~ t z is 
identified such that for length scales r « the persistent spins form a fractal with 
dimension df, for length scales r >> the distribution of persistent spins is homo- 
geneous. The zero-temperature persistence exponent, 6>, is found to satisfy the scaling 
relation 9 = Z(2 - d f ) with 9 = 0.209 ± 0.002, Z = 1/2 and d f ~ 1.58. 
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The 'persistence' problem has attracted considerable interest in recent years [1-9]. In its 
most general form, it is concerned with the fraction of space which persists in its initial 
state up to some later time. 

Hence, in the non-equilibrium dynamics of spin systems at zero-temperature we are inter- 
ested in the fraction of spins, P(t), that persist in the same state as at t = up to some 
later time t. For the pure ferromagnetic two-dimensional Ising model, P(t) has been found 
to decay algebraically [1-4] 

P(t) ~ t~ e (1) 

where 6 = 0.209 ± 0.002 [5]. Similar algebraic decay has been found in numerous other 
systems displaying persistence [9]. Most of the recent theoretical effort has gone into 
obtaining the numerical value of 9 for different models. 

Very recently, Manoj and Ray [10] have studied the spatial correlation of persistent sites 
in the Id A + A ^ model. They found that the set of persistent sites in their Id model 
forms a fractal over sufficiently small length scales. 

In this letter we present the results of an extensive numerical study of the spatial distri- 
bution of persistent spins in the pure 2d Ising model at zero-temperature. As we will see, 
the 2d Ising model exhibits behaviour very similar to that found by Manoj and Ray [10] 
in their simple Id model. 
The Hamiltonian for our model is given by 

H = -J2 SiSj (2) 

<ij> 

where Si = ±1 are Ising spins situated on every site of a square lattice with periodic 
boundary conditions; the summation in Eqn. (2) runs over all nearest-neighbour pairs 
only. 

The data presented in this work were obtained for a lattice with dimensions 1000 x 1000 (= 
AO- 
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Each simulation run begins at t = with a random (±1) starting configuration of the spins 
and then we update the lattice via single spin flip zero-temperature Glauber dynamics [5] . 
The rule we use is: always flip if the energy change is negative, never flip if the energy 
change is positive and flip at random if the energy change is zero. 
For each spin Si we define 

n i (t) = (S i (t)S i (0) + l)/2. (3) 

Hence, if rii(t) = 1 for all £ > spin Si is persistent at time t; rii(t) = otherwise. 

The total number, n(t), of spins which have never flipped until time t is then given by 

n(t) = J2i n i(t)i an d the persistence probability by [1] 

P(t) = J2<n l (t)> /N (4) 

i 

where < . . . > indicates averages over different initial conditions and histories. We averaged 
over at least 100 different initial conditions and histories for each run. 
To investigate the spatial correlations in this model, we follow Manoj and Ray [10] and 
study the 2-point correlator defined by 

C(r,t) =< ni(t)n i+r (t) > / < m(t) >, (5) 

where < . . . > now also includes the average over the lattice shown explicitly in Eqn (4). 
C(r, t) is simply the probability that spin rii+ r (t) is persistent given that rii(t) is persistent, 
averaged over the entire lattice. According to [10], the 2-point correlator satisfies the 
following dynamic scaling relation 

C(r,t) = P(t)f(r/Z(t)) (6) 

where £(£) is the persistence correlation length and f(x) is a scaling function such that 

ft \ / x~ a , for x « 1; ,„s 
f{x) ~ 1 1, for x » 1. (7) 
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As a consequence, the expected behaviour of C(r, t) in the two limits is given by 

C ^~\t-0 for r »£(*). (8) 

Clearly, as P(t) ~ t~ e , we must also have £~ a ~ t _e to satisfy Eqn (8) in the limit 
r << £(£). Assuming a power-law divergence for the persistence correlation length with t 
i.e. £(£) ~ t z then leads to the scaling relation Za = 9. As we are working with the pure 
2d Ising model at zero-temperature, we expect [11] Z = 1/2; our results are completely 
consistent with this assumption. 

To examine the correlated region (r << £(£)) we study the average number of persistent 
spins, n(/,t), in a square grid with dimensions I x I. As 

n(l,t) = [ C(r,t)rdr (9) 
Jo 

we have that 

n ,s fl 2 ~ a for Z «£(£); , . 

n(M) ~ 1 / 2 P(t) far i » (10) 

Hence, we expect the persistent spins to form a fractal with dimension df = 2 — a for 
length scales / << £(£); the distribution is homogeneous on longer length scales, namely 
for / >> £(£). We expect the crossover to occur at I ~ £(£) ~ t 1 / 2 . The scaling form for 
n(7, t) is given by 

n(l,t) = l 2 P(t)g(l/m, (11) 
where (7(x) is a scaling function satisfying 

g(x)~{ x ~ a forx«l; (12) 
V ; 1 1 for x » 1. V 7 

We now discuss our results. 

Figure 1 shows a plot of the scaling function f(x)(= C(r,t)/P(t)) against x = r/^{t) for 
various different values of t. We have assumed that ~ t 1 / 2 . The data in Fig 1 ranges 
over almost three orders of magnitude and is clearly consistent with this assumption. The 
large x behaviour of f(x) clearly follows the expected behaviour given in Eqn (7). 
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To extract a value for a we re-plot the data shown in Fig 1 on a log- log scale in Fig 2. 
The algebraic behaviour for x « 1 of the scaling function is confirmed by the linear fit. 
The slope of the straight line implies a value of a = 0.428 ± 0.007. Hence, the scaling 
relation would suggest that 9 = Za = 0.214 ± 0.004. This is, of course, consistent with 
value (0.209 ± 0.002) quoted above for 9 [5]. 

We investigate the correlated regions by obtaining a direct estimate of the fractal dimension 
df. This is undertaken by first partitioning the lattice into square grids of size / x / with I 
ranging from 4 to 250. The average number of persistent spins in each / x / square is then 
obtained. 

In Figure 3 we plot \nn(l,t) versus In I for t = 10 2 , 10 3 , 5 x 10 3 and 10 4 . We notice that 
for each of the values of t, the behaviour over sufficiently small (typically, / << yji) length 
scales is consistent with a fractal dimension df = 2 — a ~ 1.58; over longer length scales 
(typically, / >> y/i) we retrieve homogeneous behaviour (df = d = 2). Actual values of df 
range from df(t = 10 2 ) ~ 1.62 to df(t = 10 4 ) ~ 1.58. The straight lines, with slopes 1.58 
and 2.00, shown in Fig 3 are linear fits to the behaviour in the two respective regimes for 
t = 10 4 . 

We obtain an independent estimate for the exponent a by re-plotting the data for t = 
5 x 10 3 and 10 4 in scaling form. Figure 4 shows a log-log plot of the scaling function 
g(x) = n(l 7 t)/l 2 P(t) against x where x = l/y/t. We see that the data clearly fall onto a 
single scaling curve consistent with the expected behaviour given in Eqn (12). On fitting 
all of the data for \nx < —0.5 we get a value of a ~ 0.438. However, restricting the linear 
fit to Inx < —1, as indicated by the straight line in Fig 4, would imply a value of a ~ 0.50. 
Although this is slightly higher than the value we obtained from the analysis of the scaling 
behaviour of the 2-point correlator (see Eqn (8)), it is, nevertheless, consistent with our 
value of the fractal dimension in the correlated regime. 

To conclude, we have investigated the spatial distribution of persistent spins at zero- 
temperature in the pure two-dimensional Ising model. We find that the persistent spins 
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form a fractal with dimension df ~ 1.58 for length scales r << where ~ t z is the 
persistence correlation length. Furthermore, the persistence exponent satisfies the scaling 
relation 9 = Z(2- d f ) with Z = 1/2. 
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FIGURE CAPTIONS 



Fig. 1 

A plot of the scaling function f(x)(= C(r, t)/P(t)) against x where x = rfyft for t ranging 
over approximately three orders of magnitude. 

Fig. 2 

A re-plot of the data shown in Figure 1 on a log-log scale. The straight line implies a value 
of a = 0.428 ±0.007. 

Fig. 3 

A log-log plot of n(l,t) against /. Here, n(/,t) is the average number of persistent spins 
in a square (/ x /) grid at time t. The data is shown for (top) t — 10 2 ^, 10 3 +, 5 x 10 3 □ 
and 10 4 x (bottom). There is a clear crossover at / « \ft from a fractal distribution with 
dimension df ~ 1.58 to a homogeneous one with df = d = 2. The two straight lines (with 
slopes 1.58 and 2.00) are fits of the data in the two extreme cases for t = 10 4 . 

Fig. 4 

A plot of lng(x) against Inx = ln//£(t). Here the scaling function g(x) = n(/,t) /l 2 P(t). 
The straight line has slope = —0.50 and implies a value of a ~ 0.50. 
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